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ABSTRACT 

Lyman-a emitting galaxies can be used to study cosmological reionization, 
because a neutral intergalactic medium scatters Lyman-a photons into diffuse 
halos whose surface brightness falls below typical survey detection limits. Here 
we present the Lyman-a emitting galaxy LAE J095950. 99+021219.1, identified at 
redshift z = 6.944 in the COSMOS field using narrowband imaging and followup 
spectroscopy with the IMACS instrument on the Magellan I Baade telescope. 
With a single object spectroscopically confirmed so far, our survey remains con- 
sistent with a wide range of IGM neutral fraction at z 7, but further obser- 
vations are planned and will help clarify the situation. Meantime, the object we 
present here is only the third Lyman-a-selected galaxy to be spectroscopically 
confirmed at z > 7, and is ~ 2-3 times fainter than the previously confirmed 
z ~ 7 Lyman-a galaxies. 

Subject headings: galaxies: high-redshift — dark ages, reionization, first stars 

1. Introduction 

Lyman-a emitting galaxies provide a valuable probe of reionization, because resonant 
scattering of Lyman-a photons in the intergalactic medium (IGM) can suppress the observed 
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Lyman-a line by a factor of > 3 for any neutral fraction > 50% (Malhotra & Rhoads 2004 



Santos 2004). Such strong flux suppression will cause a change in the Lyman-a luminosity 

especially since Lyman-a galaxies show little evolution 



function that should be obvious- 
a t 3 < z < 6, either in luminosity function (jDawson et al. 20071 Zheng et al. 20121) or 



physical properties (Malhotra et al. 2012). Early studies concluded that the IGM neutral 



fraction was already small at z ~ 6.5 (Malhotra & Rhoads 2004 Stern et al. 2005). More 



recent work has established the Lyman-a luminosity function at both z ~ 5.7 and 6.5 to 



considerable accuracy (Ouchi et al. 2010; Hu et al. 2010; Kashikawa et al. 2011), showing 



a modest but statistically significant difference between these observed Lyman-a luminosity 
functions: The z = 6.5 luminosity function is below that at z — 5.7, and the difference can 
be adequately characterized by a pure luminosity evolution by a factor of ~ 1.3 (lOuchi et al. 



2010). 



Yet, other Lyman-a based tests for reionization — including the apparent spatial clus- 



tering of Lyman-a galaxies (McQuinn et al. 2007), the minimum ionized volume around 



observed Lyman-a sources (Malhotra & Rhoads 2006), and Lyman-a line profiles (Hu et al. 
2010 Ouchi et al. 2010) — show little evidence for neutral gas at z ~ 6.5. This leaves an 
open question — is lower Lyman-a luminosity function at z = 6.5 due to neutral gas, or is it 
an intrinsic evolution in the galaxy populations? 

To distinguish between these possibilities, we can look to still higher redshifts, where 
the IGM neutral fraction should be higher and its effects on Lyman-a stronger. The highest 
redshift readily accessible to Lyman-a searches using CCDs is z ~ 7.0, in the 9650A window 
in the night sky OH emission spectrum. We are pursuing a 9650A narrowband survey using 
the IMACS imaging spectrograph on the 6.5m Magellan I Baade Telescope at Las Campanas 



Observatory (Hibon et al. 2011). We surveyed 465 square arcminutes, corresponding to 
~ 72000Mpc 3 . After a careful selection, we found 6 z~6.96 LAE candidates (Hibon et al. 



2011). To confirm whether these are real LAEs, we obtained multi-object spectra with 



IMACS. In this Letter, we present the spectrum of LAE J095950. 99+021219.1, which was 
identified as a candidate redshift z ~ 7 Lyman-a emitting galaxy (candidate LAE 3) in 
Hibon et al. (2011). Our spectroscopy reveals a single, isolated Lyman-a line at redshift 



z = 6.944. 

Throughout the paper, we adopt a ACDM "concordance cosmology" with Qm 
Q A = 0.73, and H = 71 kms" 1 Mpc" 1 . 



0.27, 
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2. Spectroscopic Observations and Analysis 
2.1. Observations 

We observed our candidate z ~ 7 Lyman-o; galaxies using the Inamori Magellan Areal 
Camera and Spectrograph (IMACS) on the 6.5m Magellan I Baade Telescope on the nights 
of 29-30 December 2010, and 8 February 2011. The February data were of lower quality and 
are not used here. We used custom multi-slit masks, shared between two primary observing 
programs. We selected the f/2 camera and the 300-line red-blazed grism with 1" slitlets as 
the best compromise between areal coverage, spectral coverage, and spectral resolution. 

Observations were split among five slit masks (two per night in December, and one 
in February). The time per mask and observing conditions are summarized in table [I] 
While the position angle of the masks were not all identical, the data were taken with- 
out dithering the telescope. Moreover, the targets were centered on their slitlets, and 
LAE J095950. 99+021219.1 is compact compared to the seeing. This allows us to combine 
all of the spectroscopic data into a single ID spectrum (see below). 



2.2. Data Reduction 

We performed initial data reduction steps using the COSMOS software package^ COS- 
MOS steps include bias frame subtraction, spectroscopic flat fielding using continuum (quartz) 
lamp exposures, and wavelength calibration using arc lamp exposures. COSMOS also sky- 



subtracts the spectra, using the Kelson (2003) algorithm to remove night sky lines. Finally, 



COSMOS extracts a 2D, rectified spectrum for each slitlet. 



2 Note the double meaning of the acronym COSMOS. We deny any responsibility for the ensuing confusion. 



Mask ID 

COSMOS1 
COSMOS2 
COSMOS3 
COSMOS4 
COSMOS-Feb 8 Feb 2011 



Observation 
date (UT) 
30 Dec 2010 

30 Dec 2010 

31 Dec 2010 
31 Dec 2010 



Number of 
exposures 

3 
4 
3 



Time per 

exposure 

1800 

1800 

1800 

1860 

1800 



Seeing 

(approx) 

0.5-0.9" 

0.5-0.9" 

1-1.5" 

1-1.5" 

1.2-2" 



Table 1: Summary of slit mask observations. 
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We performed subsequent steps in two ways, either (a) combining exposures from each 
mask separately, and then combining results from different masks; or (b) directly combining 
all exposures from multiple masks. 

Treating masks separately gives four 2D spectra of LAE J095950. 99+021219.1 from 
December, and one more from February. Most of these 2D stacks show a weak but visible 
emission line in the spectrum of LAE J095950. 99+021219.1. We next combined the four 
December spectra into a stacked 2D spectrum comprising our best 7.05 hours of data. To do 
this, we first averaged the four 2D spectra. Next, we made a median-combined stack. We then 
subtracted the two, and computed the (sigma-clipped) noise level in the difference. Finally, 
we constructed a hybrid stack, using the value from the average stack almost everywhere, but 
the value from the median stack wherever the difference between these two stacks exceeded 
10cr. This yields a lower noise estimate than the median, yet remains more robust to outliers 
than the mean. The emission line becomes readily evident in this combined stack. 

To test the robustness of our results, we also combined all December exposures in single 
14-frame stacks, using various outlier rejection schemes (median stacking with 3- and 5cr 
rejection, and average stacking with 2.5a rejection). The emission line remains comparably 
significant in all of these stacks. The stacked 2D spectrum around the emission line, using 
average stacking and 2.5a rejection, is shown in figure [1} 

We also made stacks by bootstrap resampling, stacking 14 exposures selected randomly 
with repetitions permitted. We remeasured the flux at the location of the detected line, 
using aperture photometry in the 2D stacks. The bootstrap fluxes were 104% ± 10% of 
the "normal" stack flux for a 10-pixel diameter (2" x 20A) aperture, and 116% ± 21% for 
a 14-pixel diameter (2.8" x 28A) aperture. Among 1000 bootstrap simulations, the lowest 
measured fluxes were 77% and 66% of "normal" for the 10- and 14-pixel apertures. The 
observed line is therefore not a fluke caused by a handful of exposures. 

We next extracted a ID spectrum from the two-stage stacking (method "a"), using the 
IRAF task "apall" with an unweighted extraction of 1.4" (7 pixel) window width, centered on 
the emission line and parallel to the dispersion axis. (COSMOS 2D output has the dispersion 
axis parallel to the x-axis, so we need not fit a trace to the continuum, which is undetected in 
the present data anyway.) We performed no further sky subtraction, since that too is done by 
the COSMOS package. To get another estimate of the noise level in the data, we extracted 
five further ID spectra from the 2D stacked spectrum, each at a different spatial position 
along the slit. Each should be essentially a pure noise spectrum. The variance among these 
five parallel traces provides a wavelength-dependent noise estimate, and the line is significant 
at the 4.5a level against this estimate. The extracted ID spectrum is shown in figure [2j 



-5 - 



Significance: To explore the significance of the line detection, we measured aperture fluxes 
at a grid of clean locations in the 2D spectrum of LAE J095950. 99+021219.1 (after rescaling 
the 2D spectrum by the noise ratio cr(9658A)/<r(A)). The RMS counts among these apertures 
corresponds to a la noise of 1.34 x 10~ 18 erg cm -2 s -1 , against which our line is a 6.3<r event. 
Among > 400 non-overlapping 5-pixel apertures, the brightest two were 4.1er and 3.2<r events 
(65% and 51% as bright as the LAE J095950. 99+021219.1 line), suggesting only mildly non- 
gaussian noise. 

The search that found this line was based on 6 candidates, each with a position known 
to 1" and an expected line wavelength known to AA = 90A. Given our spatial and spectral 
resolution, this corresponds to ~ 6 x 2 x 20 = 240 independent resolution elements. Our 
significance level estimates range from 4.5a to 6.3<r, and for Gaussian noise, the corresponding 
chance probabilities in 240 trials range from < 10~ 3 to < 10~ 7 . 



Spectroscopic flux calibration: Each mask included two blue stars with well measured 



photometry from the COSMOS project (Capak et al. 2007). We flux calibrate the observed 



emission line of LAE J095950. 99+021219.1 by direct comparison with the observed counts in 
one of these stars, which was observed under identical conditions as our science targets. Both 
the emission line count rate and the comparison star's count rate per unit wavelength were 
measured directly in the 2D spectra and at the same wavelength. For the emission line, we 
used the 10-pixel diameter (2" x 20A) aperture. This yields a flux of 8.5 x 10~ 18 ergcm~ 2 s -1 
for LAE J095950. 99+021219.1, corresponding to line luminosity 4.9 x 10 42 ergs _1 at z — 
6.944. 

The fractional uncertainties in this flux are ~ 20% from photon counting statistics (with 
a statistical signal-to-noise ratio of s/n = 4.5-6.3 from the ID or 2D spectra), 12% from the 
choice of aperture used to measure the line flux in the 2D spectrum, and < 10% from the 
assumption that the comparison star's flux density at 9680A equals its z-band flux density. 
The final spectroscopic flux measurement is (8.5 ± 2) x 10~ 18 erg cm -2 s . 



2.3. Comparison to narrowband imaging results: 



The observed spectroscopic line flux is smaller than the narrowband flux (Hibon et al. 



2011). Part of the discrepancy can be attributed to continuum in the narrowband filter. The 



emission line is near the blue edge of the filter bandpass, so continuum emission redward of 
the line will contribute relatively strongly to the narrow band flux. 



We also re-examined the narrowband flux measurements for LAE J095950. 99+021219.1. 
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Following Hibon et al. (2011), we used moderately bright stars that are well detected but 



unsaturated in both the public COSMOS z-band image and our NB9680 image. The nar- 



rowband magnitudes in Hibon et al. (2011) were based on 1" diameter aperture fluxes, with 



an aperture correction based on the difference between SExtractor "magiso" and aperture 



fluxes (see Hibon et al. (2011) for more details). In the present work, we omitted the aper- 
ture correction step, using instead identical 1.2" diameter flux measurements for both the 
science objects and the reference stars in the NB9680 image. (The precise aperture di- 
ameter is unimportant, since LAE J095950. 99+021219.1 is compact compared to the point 
spread function.) We obtained a narrowband magnitude N B96S0 ab = 24.86 ± 0.18. The 
corresponding narrowband flux is 12 x 10~ 18 ergcm~ 2 s _1 . 

This is almost 50% more than the spectroscopically determined emission line flux (a 
1.5cr difference). If we attribute the difference to continuum emission in the narrowband 
filter, the corresponding flux density is fx ~ (/nb — /s P ec)/AA ps 3.5 x 10~ 18 /60A ps 5.8 x 
10~ 20 ergcm -2 s _1 A -1 , or AB = 25.75. The implied observer-frame equivalent width is 
~ 120A (though consistent with an arbitrarily large equivalent width at 1.5c). While the 
object is undetected in the z' filter down to AB > 26.4, a continuum magnitude of 25.75 
redward of the line at 9651 A remains allowed, since most of the z' filter's transmission lies 
blueward of that wavelength. The object is also undetected in the WIRDS J-band image 



dBielby fc al|2012, in prep[ ), with a 1.2" aperture flux ~ (1.8 ±2.1) x lO^ergcm" 2 s" 1 Hz" 



corresponding to Jab ^ 24.1 (3a). 



3. Interpretation 

We interpret the line at 9657A in LAE J095950. 99+021219.1 as Lyman-a at redshift 
z = 6.944, based on non-detections in all filters blueward of this line, and on the absence of 
other optical lines. 

Were the primary line Ha (at z ps 0.472), or [O III] A5007 (at z ps 0.928), we would ex- 
pect other prominent emission lines in our IMACS spectrum. Figure [3] shows non-detections 
in the ID spectrum at the expected line wavelengths. Corresponding upper limits are 
summarized in table [2j The "Ha" case (z ~ 0.472) is disfavored by the non-detections 
of [O III] A5007 and [O II] A3727, with 3a line ratio limits /(O III) //(Ha) < 1/2 and 
/(O II)//(Ha) < 2/3. If the primary line is [O III] A5007, unfortunately placed night-sky 
line residuals overlap the expected locations of [O III] A4959 and H/3, precluding interesting 
limits. Fortunately, the expected [O II] A3727 line location is clean, and gives a tight upper 
limit /(O II)//(0 III) < 0.2 (3a). This provides some evidence against the [O III] A5007 
interpretation. Ratios of /(O II)//(0 III) this small are seen in a significant minority of star 
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forming galaxies (Xu et al. 2007 McLinden et al. 2011 Richardson et al. 2011, in prep; Xia 



et al. 2012, in prep), but more can be ruled out by this line ratio. 



To address the [O II] A3727 possibility and further improve our constraints on [O III] 



A5007, we examine equivalent widths. Following Hibon et al. (2011), we combine our spec- 
troscopic line flux with optical magnitude limits of 27.9, 27.6, and 27.3 mag (5a, AB) in 
g', r', and i' filters respectively. Since star-forming galaxies have f v ~ constant, we have 
/ A (9657A) w 3.6 x lO" 20 x i(r a4x28 x c 

detection in these optical images then implies a 5a limit EW = f ]in J fx 
10~ 21 )A = HOOA (observer frame). 



/(9657A) 2 = 7.3 x 10" 21 ergcm- 2 s~ 1 A" 1 . The non- 

(8 x 10~ 18 /7.3 x 



While [O III] A5007 and Ha emission line sources with equivalent widths this large exist 



(Rhoads et al. 2000; Kakazu et al. 2007; Straughn et al. 2008, 2009; van der Wei et al. 2011 



Atek et al. 2011), they are exceptional, rare objects. In Hibon et al. (2011), we estimated 
the numbers expected in our survey based on published line luminosity functions (Kakazu 
et al.||200"7 Geach et al.||2010 ) and equivalent width distributions ( Straughn et al.|[2009 ). We 
found that < 0.3 [O III] A5007 emitters and < 0.6 Ha emitters are expected. 

[O II] A3727 emitters have generally smaller equivalent widths. We found no [O II] 



A3727-selected objects with EW > 1100A/(1 + Z[oir\) ~ 425A in the samples from Straughn 



et al.| Q2009D (30 objects), |Kakazu et aL| Q2007D (24 objects), |Xia et al.| ( |2012, in prepj ) (11 
objects), or prozdovsky et al.| ( pOOlj ) (400 objects). Thus, < 1/465 = 0.0022 of [O II] 
A3727 emitters might enter our sample as LAE candidates. The luminosity function from 
Rigopoulou et al. (2005) suggests that our survey volume should contain ~ 45 [O II] A3727 



galaxies. Among these, < 0.1 object should pass our Lyman-a selection criteria. 

Thus, the aggregate sample of foreground emitters expected in our survey is < 1 galaxy. 
In contrast, our survey volume at z ~ 6.95 should contain between ~ 2.5 and ~ 11 Lyman-a 
galaxies with line fluxes > 8 x 10 -18 erg cm -2 s _1 , based on the z ~ 6.5 luminosity functions of 



Hu et al. 


(2010 


) and 


Ouchi et al. 


(2010 



of the observed emission line. 



4. Discussion 

Redshift z ~ 7 is the current frontier in reionization studies, an area of active exploration 
where our observational knowledge is growing rapidly. 

The recently discovered quasar at z ~ 7.1 (Mortlock et al 2011), combined with spectro- 
scopic followup (and occasional confirmation) of z ~ 7 galaxy candidates from HST WFC3 
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surveys, provide an unprecedented look at this epoch. Observations of these objects seem to 
favor the continued existence of significant neutral intergalactic gas as late as z ~ 7. This 
is surprising, given that microwave background polarization data from WMAP favor a char- 
acteristic reionization redshift z r ~ 11, and that the IGM at z ~ 6.2 is highly ionized, with 
a neutral fraction of only 1-4% based on quasar spectra (Fan et al 2006). Nonetheless, the 
ionized bubble around the z ~ 7.1 quasar appears too small to be comfortably explained in 
a fully ionized medium, unless the quasar is itself remarkably young (~ 10 6 years) (Bolton et 
al 2011). Similarly, three independent research groups have argued that the fraction of z ~ 7 
Lyman break candidates showing Lyman-a emission appears smaller than would be expected 
in an ionized medium (Vanzella et al 2010; Stark et al 2010; Ono et al 2011; Schenker et 
al 2011; Pentericci et al 2011). Still, these results depend on the reliability of photometric 
selection criteria, and a contamination of order 50% could explain the observations without 
recourse to neutral gas (e.g., Schenker et al 2011). 

Lyman-a galaxy surveys offer a complementary approach to studying reionization. The 
underlying physics is the same as for spectroscopic followup of Lyman break samples, but 
the survey selection proceeds differently, leading to different potential selection biases. The 
uncertainties in the method are likewise very different from those associated with the quasar 



near-zone measurement (Mortlock et al. 2011 Bolton et al. 2011), or the Gunn-Peterson 



trough (Fan et al. 2006). Because of this, conclusions about cosmological reionization will 



be strongest when they are based on multiple independent methods. 

The work we present here is only the second large-area narrowband survey for for Lyman- 
a galaxies at redshift z m 7.0, following on the work of lye et al. (2006) and Ota et al. (2010). 



The spectroscopic confirmation of LAE J095950. 99+021219.1 demonstrates that such objects 
can be identified at flux levels considerably fainter than the 2 x 10~ 17 erg cm~ 2 s _1 line of IOK- 



1 (lye et al. 2006) or the recently reported z ~ 7.2 narrowband-selected Lyman-a galaxy 



SXDF-NB1006-2 (Shibuya et al. 2011). Our observed line flux, 8.5 x 10" 18 erg cm" 2 s -1 , 



corresponds to a rest-frame line luminosity of only Li ya = (4.9 ± 1) x 10 42 ergs _1 . This is 
at or below the characteristic line luminosity L* Lya from Schechter function fits to z ~ 6.5 
Lyman-a samples (e.g., L* Lya = (4.4 ±0.6) x 10 42 ergs _1 and $* = 8.5 x 10~ 4 cMpc~ 3 , 
(|2010[); L* Lya = 1.0 x lO^ergs" 1 and $* = 6 x IQ~ 5 cMpc~ 3 



et al. 



Hu et al. 



(2010)). 



Ouchi 



Thus, our survey has achieved sensitivity to typical Lyman-a emitters. Were the z = 7 
luminosity function unchanged from that at z = 6.5, we would expect our survey volume 
to contain 11 or 2.6 Lyman-a emitters brighter than LAE J095950. 99+021219.1, based on 



the LF of Ouchi et al. (2010) or Hu et al. (2010) respectively. This reflects large differences 



in expectations from different published z = 6.5 luminosity functions. We will present a 
detailed analysis of our survey's constraints on the Lyman-a luminosity function in a future 
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paper, after we have more sensitive spectra for our remaining candidates. For now, the 
galaxy LAE J095950. 99+021219.1 is among the few most distant spectroscopically confirmed 
galaxies known, and the faintest to be discovered through a direct search for Lyman-a line 
emission. 
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Fig. 1.— 2D spectra of LAE J095950.99+021219.1. Upper panel: Full December data set 
(14 exposures, 7.05 hours' integration) stacked together. Wavelength increases from 9550A 
at left to 9800A at right. The spatial extent is 12" from bottom to top. Lower panels: First 
night (left) and second night (right) of data separately. The lower panels have been lightly 
smoothed for clarity, with a Gaussian kernel having a = 0.6 pixel. Cyan bars at top and 
bottom mark the wavelengths of night sky emission lines, with longer bars denoting brighter 
lines. 
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Fig. 2. — The extracted ID spectrum (lightly smoothed with a "1 2 1" filter, with 4A 
FWHM) at top, and a noise estimate (unsmoothed) at bottom. The sole significant feature 
is the Lyman-a line at 9657A, corresponding to redshift z = 6.944. Some residuals from 
subtraction of the strong night sky OH emission lines remain, and are marked. The most 
notable of these is the spike at 9698A. The plotted noise is based on the sigma image from 
the 14-frame stack, scaled for the number of spatial pixels combined in the ID extraction. 
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Formal flux 


3<r limit 
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(10 18 erg cm 2 s x ) 
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[0 III] A4959 


7299A 


-0.6 ±0.8 


1.8 


0.472 


H/3 


7154A 


-0.2 ±0.8 


2.2 


0.472 


[0 II] A3727 


5486A 


0.6 ± 1.6 


5.4 


0.928 


[0 III] A4959 


9565A 


-4.3 ±4.6 


9.5 


0.928 


H/3 


9375A 


—12 ± 15 


33 


0.928 


[0 II] A3727 


7188A 


0.0 ±0.5 


1.5 



Table 2: Line flux limits at locations where we would expect emission lines, were the 
primary line at 9657A not Lyman-a. z = 0.472 corresponds to Ha at 9657A, and 
z = 0.928 corresponds to [O III] A5007 at 9657A. For comparison, the 9657A line has flux 
8.5 x 10~ 18 ergcm~ 2 s _1 . Limits tabulated here are based on fitting a Gaussian line profile 
to the observed spectrum, with the central wavelength and line width for the fit fixed by the 
properties of observed line. 
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5450 5500 5550 7200 7300 

Wavelength (A) Wavelength (A) 



Fig. 3. — If the emission line in LAE J095950. 99+021219.1 were not Lyman-a but instead 
were Ha or [O III] A5007, we should expect additional emission lines in our optical spectrum. 
These figures show the expected locations of the [O III] AA4959,5007, H/3, and [O II] A3727 
lines, for the cases where the primary line is either [O III] A5007 at z — 0.929 (top two 
panels) or Ha at z — 0.472 (bottom two panels). In no case do we see a prominent emission 
line, although for the z = 0.929 case the [O III] A4959 and H/3 lines both fall on night sky 
line residuals, leaving the absence of an [O II] A3727 line at 7188A and the lack of any optical 
broad band detection as the best discriminants against this possibility. As in fig. 2, the signal 
spectra are lightly smoothed with a (1 2 1) filter, while the noise spectra are unsmoothed. 



